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We propose two sets of initial conditions for magnetohydrodynamics (MHD) in which both the 
velocity and the magnetic fields have spatial symmetries that are preserved by the dynamical equa- 
tions as the system evolves. When implemented numerically they allow for substantial savings in 
CPU time and memory storage requirements for a given resolved scale separation. Basic properties 
of these Taylor-Green flows generalized to MHD are given, and the ideal non-dissipative case is 
studied up to the equivalent of 2048'^ grid points for one of these flows. The temporal evolution of 
the logarithmic decrements 5 of the energy spectrum remains exponential at the highest spatial res- 
olution considered, for which an acceleration is observed briefly before the grid resolution is reached. 
Up to the end of the exponential decay of 5, the behavior is consistent with a regular flow with no 
appearance of a singularity. The subsequent short acceleration in the formation of small magnetic 
scales can be associated with a near collision of two current sheets driven together by magnetic 
pressure. It leads to strong gradients with a fast rotation of the direction of the magnetic held, a 
feature also observed in the solar wind. 
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I. INTRODUCTION 

Astrophysical and geophysical flows arc highly turbu- 
lent, with strong coupling between a wide range of spatial 
and temporal scales. The complex behavior of such flows 
is far from understood, and their study through direct 
numerical simulation (DNS) in three space dimensions is 
limited to modest scale separation, even at the largest 
resolution achieved today [H, 0]. In fact, many numerical 
experiments designed to study turbulent flows are con- 
fined to Cartesian geometry and periodic boundary con- 
ditions in order to maximize the scale separation (or the 
Reynolds number in the dissipative case) in the compu- 
tations. Furthermore, scaling laws of turbulence param- 
eters, such as the incrtial index of the energy spectrum, 
the Kolmogorov constant, or the skcwncss and flatness 
of velocity derivatives, all converge slowly to their high 
Reynolds number limits (assuming they exist), and they 
converge at different speeds Q. 

There are several ways to obtain a larger scale separa- 
tion in simulations of hydrodynamics and magnetohydro- 
dynamics (MHD). Adaptive mesh refinement (AMR) is 
a method of improving resolution locally, based on struc- 
tures that develop in the flow, without having to refine 
the grid everywhere. For dissipative turbulence, however, 
AMR can be difficult or costly to implement, unless the 
emerging structures are strongly concentrated in space 
Q . Success with AMR has been achieved in certain con- 
texts, such as in the study of the possible development 
of singularities in hydrodynamics and MHD 0, Q . 

In more homogeneous flows, an alternative is to model 



the small scales that are not explicitly resolved in a simu- 
lation. Such methods abound for neutral fluids (see e.g., 
7 1), but in the case of MHD the research is less developed 
8| . An additional difficulty in MHD is the lack of labo- 
ratory experiments at large magnetic Reynolds number 
to evaluate such models against DNS and thus to fur- 
ther improve them Also, although they are useful 
for many applications, these methods cannot be used to 
study mathematical problems such as the development 
of singularities in finite time, since the original equations 
are, as a rule, modified or replaced by a different set of 
equations to mimic the effect of the turbulent mixing and 
dissipation at unresolved scales. 

A third approach is to construct a flow with symme- 
tries that are hypothesized to hold under the evolution 
of the relevant dynamical equations. Numerical simu- 
lations of such flows have been performed using codes 
designed to exploit their symmetries, leading to signifi- 
cant advances in the study of turbulent flows [l^, El , due 
to substantial savings in both computing time and mem- 
ory (see also [HI). Like AMR, flows with symmetries are 
suitable for studying t he p ossible development of singu- 
larities in an ideal flow [ij] . Efficient codes implementing 
fiow symmetries have also been utilized in the study of 
MHD flows, although in the context of a forced turbu- 
lent dynamo rather than singularity [3, [IB| • It is worth 
noting that in both of these cases, although the magnetic 
field was evolved by a code employing symmetries, it was 
initialized as a small, random seed field whose symmetries 
were incidental and not deterministic; therefore, only a 
symmetric velocity field needed be prescribed. 
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In this paper we propose two paradigmatic MHD flows 
with symmetries. Numerical implementation of these 
symmetries are used to demonstrate that the symme- 
tries are preserved by the MHD equations as the system 
evolves for many turnover times, as evinced by a com- 
parison of simulations with and without imposed symme- 
tries. As the implemented symmetries allow for substan- 
tial savings in memory and computing time, we are able 
to study the evolution of an ideal MHD flow at unprece- 
dented spatial resolution. The structure of the paper is 
as follows. In the next section we present the equations 
and definitions of the possible flows that can be studied 
in this context, giving as well some of their basic geo- 
metrical properties. Section |lll] is devoted to a study of 
the ideal non-dissipative case and to the rapidity of the 
development of small scales through an analysis of the 
temporal evolution of the analyticity strip (to be defined 
below). At the end of Section Hill we also discuss numer- 
ical simulations showing that the symmetries of the fiow 
are preserved for all times considered in the present work. 
Finally, Section [TVl presents our conclusions. 

II. TAYLOR-GREEN FLOWS 

A. The MHD equations 

We recall here the ideal MHD equations for an incom- 
pressible flow; these equations stem from the conserva- 
tion of momentum and from Maxwell's equations for ve- 
locities substantially less than the speed of light. In terms 
of a dimensionlcss velocity v and magnetic induction b, 
they are 

dtv + v-Vv = -VP+jxb, (1) 
dth = V X (v X b) , (2) 

together with V-v = = V-b;j = Vxbis the 
current density and P the pressure. Note that b is in 
fact expressed in terms of an Alfven velocity. The total 
energy 

E=^{w^ + h^)=Ev+EM, (3) 

(where Ey is the kinetic energy, Em is the magnetic en- 
ergy, and ( ) signifies integration over the domain), the 
cross-correlation 

^c = ^(v-b), (4) 

and the magnetic helicity 

HM = l{a-h), (5) 

(where b = V x a, and a is the vector potential) are 
conserved by these equations in three space dimensions. 
The conservation of these quantities (up to the numeri- 
cal truncation error) is monitored in the simulations and 
used as one criterion to decide the last reliable time for 
which the accuracy of the solutions is satisfactory. 



B. Hydrodynamic Taylor— Green flows 

The Taylor-Green (hereafter, TG) vortex is an im- 
portant hydrodynamic flow for the insights it provides 
at moderate complexity but low cost. It mimics von 
Karman flows between two counter-rotating disks, as 
used in several turbulence experiments (including recent 
experiments to reproduce generation of magnetic fields 
by dynamo action; see e.g.[l6j and references therein). It 
was originally motivated as an initial condition that, in 
spite of its symmetries, would lead to the rapid develop- 
ment of small spatial scales and therefore proposed 
as a paradigm of the direct cascade of energy in a turbu- 
lent fiow. Analytic expansions were considered to s tudy 
the possible development of singularities in this fiow [3 ■ 
Fundamental to the TG vortex in later numerical stud- 
ies were its symmetry properties that allowed for gains 
of a factor of four in linear resolution [l^, [l^, [2l| , a 
feature that was also exploited to study the generation 
of magnetic fields by dynamo action [lj| . 

The simplest TG vortex is given by the velocity field 
initial data 

Vx = Wo sin a; cos y cos z (6) 
Vy = — Wo cos X sin 2/ cos z (7) 
V, = 0. (8) 

The z-component of the velocity grows with time as the 
system evolves, due to the pressure gradient in the mo- 
mentum equation. When the hydrodynamic equations 
arc integrated in time starting from this initial condi- 
tion, the symmetries of the fiow are preserved for many 
turnover times. 

Within a periodic cube of length 2tt, there are two 
planes of mirror symmetries (or anti-symmetries) in each 
dimension: x = 0, tt; y = 0, tt; and z = 0, tt. The 
fiow is also invariant by rotation of tt around the two 
axis X = z = tt/2 and y = z = 7r/2 and by rota- 
tion of 7r/2 around the third axis x — y — t:/2. As a 
consequence of the last rotational invariance, note that 
Vy{x,y,z) = —Vx{tt — y,x,z). Exploiting these symme- 
tries allows one to recover the dynamics of the entire 
domain from computations in a "symmetry" box of size 
[0,7r/2]'^ [3, [l^. By reducing the amount of informa- 
tion needed to represent the full cube, the symmetries 
allow for simulations at much higher Reynolds numbers, 
or simulations of the ideal inviscid fiow for a longer time 
than a full DNS for a given fixed cost. Specifically, at a 
given computational cost, the ratio of the largest to the 
smallest scale available to a computation with enforced 
Taylor-Green symmetries is enhanced by a factor of 4. 
Equivalently, at a given Reynolds number, a factor of 
4 in savings in linear resolution offered by the Taylor- 
Green fiow thus leads to a factor of 32 savings in total 
computational time and a factor of 64 savings in memory 
usage. 
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C. Two possibilities in MHD 

With the velocity given above, there are several possi- 
bilities for the choice of initial conditions for the magnetic 
field with the same symmetries. The first set of magnetic 
initial conditions will be called the "insulating" case, with 
b everywhere normal to the walls, and the current j con- 
tained within what can be called the "insulating box" 
[OjTt]"^; hence, the walls in this box are insulating, with 
the magnetic field behaving like the vorticity in the hy- 
drodynamic TG case: 

6^'^ = 60 cos a; sin y sin z (9) 
by ^ = &o sin X cos y sin z (10) 
b^'' ~ — 26o sinxsinj/cosz. (11) 

The second possibility in MHD is called the "conduct- 
ing" case; it has the magnetic field, just like v, contained 
entirely within the box: 

fe^"' = bo sin 2x cos 2y cos 2z , (12) 
by^"^ = foo cos 2a; sin22/ cos 2z , (13) 
b'f'^ = -26ocos2a;cos2?/sin2z . (14) 

For this flow He = (uncorrelated v and b), and 
it yields a current that is everywhere perpendicular to 
the walls of the conducting [0, tt]^ box. In this paper, we 
concentrate on the study of the first flow, which we name 
IMTG ("insulating magnetic Taylor-Green") hereafter. 

D. Basic properties of the IMTG flow 

The initial velocity in the insulating box has a shear 
layer between two counter-rotating eddies, and the mag- 
netic field is initially proportional to the vorticity. Lo- 
cal kinetic and magnetic helicities are anti- symmetric 
with respect to the planes of mirror symmetries and 
thus vanish globally. Kinetic helicity, which is defined 
as Hji = {v ■ u)) /2, is not generally conserved by the 
ideal MHD equations. 

Noting that Ey = v^/S, Em = 36g/8, He = and 
Hm = 0, in what follows we take vq = 1 and 60 = 1/V3 
(and thus Ey ~ Em = 1/8 at t = 0). Fourier space 
components of the symmetric fields are non-zero only for 
jointly even or jointly odd wavenumbers. Even and odd 
components are treated separately in the code, using spe- 
cially designed Fourier transforms that are performed on 
N/4: points, where the resolution N denotes the number 
of grid points in the full periodicity box. Computations 
are dealiased using the 2/3 rule, and the wavenumbers 
are thus integers in the range [fc,„i„ = l,kmax = -^/3]. 
Finally, time differencing is achieved using, in general, a 
second-order Runge-Kutta method. 

When the IMTG flow is used as the initial condition 
for a freely decaying simulation, the symmetries are pre- 
served by the time evolution for many turnover times. 



This is shown in the next section for the ideal case (of 
interest for the present work), and will be shown for the 
viscous case in a forthcoming work. 

III. NUMERICAL RESULTS 

A. Temporal evolution in the ideal case 

In the ideal case, E, Hm, and He are conserved (ex- 
cept for time-stepping errors), and the flow, with a finite 
number of modes, eventually evolves to statistical equi- 
libria that depend on these conserved quantities [2^ . In 
the simplest case, equipartition between modes obtains, 
leading to an accumulation of energy in the smallest 
scales with a k'^~^ spectrum in d space dimensions, the 
so-called ultra-violet energy catastrophe. It was shown 
recently [2l| that for d = 3, the small-scale fc^ spectrum 
acts as an eddy viscosity to the modes at large scales, 
thereby putting on a strong footing the assertion that 
the ideal problem contains the dynamics of the dissipa- 
tive problem, at least for times short enough that the 
equipartition has not yet been able to reach all scales. 

B. The singularity test 

In this light, it is of interest to examine the behavior of 
the IMTG flow we just proposed. One central question 
is whether ideal MHD flows develop singularities in a fi- 
nite time. In a celebrated paper [23[ it was shown that, 
for the Euler equation, the norm of vorticity con- 
trols the breakdown of solutions; this result was extended 
to MHD, in which case the singularity is controlle d by 
the combined L°° norm of vorticity and current [l^, l25l |. 
There are a number of carefully constructed numerical 
studies devoted to ideal behavior for both Euler and for 
ideal MHD, with no clear-cut conclusion of a specific flow 
showing signs of singularity (see 0, [2^, [l^ for the Euler 
case and [1, [2^ [2^ for MHD); note that the particular 
choice of initial conditions is essential, as exemplified in 

In this context, we show in Fig. [T] the evolution of the 
maximum of the current jmax in semi-log coordinates in 
a run computed on an equivalent grid of 2048'^ points. 
An approximately exponential phase up to t « 2.47 is 
followed by a more abrupt growth corresponding to a 
change in the spatial location of jmax- The exponential 
behavior at early times has been observed by a number of 
authors, and can be related to a shear-induced instability. 
In fact, an exponential fit is particularly good in the range 
t = 1.80 to t = 2.27, where the least square error ^ « 
8 X 10^'^. A power-law fit also appears appropriate in 
the range t = 2.10 to t = 2.47, with now a least-square 
error £^ ^ 2 X 10~^. The inset in Fig. [T] shows a blow-up 
of the evolution of jmax, in log-log coordinates, with this 
power-law fit up to the sharp transition at t « 2.5. In 
spite of the steepness of the index for the power-law fit 
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FIG. 1: Maximum of the current in semi-log in a 2048^ IMTG 
run; the insert gives a log-log zoom, with triangles for data, 
dash-dot for the fit with [f, — t]^^, and solid line for a power 
law. 



1 .00 










0.10 




1 78 














312 


>^ 




0.01 






■\ 

\ 












5.0 










4.5 










4.0 










3.5 










3.0 










2.5 










7.0 










0.5 


1.0 


1.5 2.0 


2,5 



FIG. 2: (Color online) Top: Logarithmic decrement of the 
total energy spectrum for various resolutions in the 128^ (long 
dash), 256^ (dash-triple dot), 512^ (dash-dot), 1024^ (dash) 
and 2048^ (solid) runs. Note the acceleration in 5 in the 2048 
run at f « 2.5 Bottom: Spectral index in the 128^ (-f), 256'^ 
(*), 512^ («.), 1024^ (A) and 2048^ (solid line) runs. 



r ), the error for the power-law fit is about four 



orders of magnitude smaller than for a jma x ~ — t]~°' 
fit in the same range, as proposed in psi . [2^ . 25j for a 
singularity to occur at a finite time t*. The inset in Fig. 
[T] also shows the best fit to this expression with a = 1 . 

The resolution of the run does not allow us to com- 
pute far beyond the time shown. Indeed, from a visual 
inspection of the flow, the computation can be consid- 
ered reliable up to t « 2.61, at which time aliasing rip- 
ples in the current appear, roughly at the 15% level (not 
shown). As a comparison, the level of aliasing in the cur- 
rent a± t K, 2.55 is 9% (below we discuss two different 



accuracy criteria based on the evolution of the energy 
spectrum and the conservation of the MHD invariants). 
We are thus led to conclude that, due to the substantially 
smaller errors of an exponential or power-law fit than for 
a [i* — fit to the norm of the current or vorticity, 
the flow shows no sign of developing singularities for the 
reliable duration of the computation. Additional studies 
will be needed to examine the later temporal evolution 
of the ideal flow at higher resolution. 



C. Logarithmic decrement and accuracy criteria 

Another way to determine whether a flow remains reg- 
ular at all times or develops a singularity is to measure 
the logarithmic decrement. This quantity is also use- 
ful as an accuracy criterion for gauging up to what time 
the ideal computation is reliable. The logarithmic decre- 
ment is defined through a (logarithmic) fit of the total 
energy spectrum E(k,t) [defined in the usual way, with 
/ E{k,t)dk = E{t)] with the expression 



E{k,t) ^ C(t)fc-"(*)e-2*(*)''' . 



(15) 



The logarithmic decrement 5 must not reach zero (e.g. 
by decaying exponentially in time, 5(t) = Soe~*^'^) for 
the fiow to remain regular at all times. Independently 
of its time dependence, the instantaneous value of S can 
also be compared with Ax (the distance between spatial 
grid points) to quantify the accuracy of the computation. 

We show in Fig. [2] the temporal evolution of 6 com- 
puted on the total energy spectrum and for simulations 
with resolutions ranging from 128'^ to 2048"^ grid points. 
A striking result that appears only in the highest-resolved 
run, again at t « 2.5, is an acceleration in the decrease 
of S. This sharp change in the evolution of S is observed 
when the fit given by Eq. is done using the total 
energy spectrum or the magnetic energy spectrum from 
the simulation, but not when using its kinetic counter- 
part. These results were verified to be independent of 
the timestep used (which was halved twice) , as well as of 
the order of the temporal scheme (which was tested at 
second and fourth order). 

Interestingly, the accelerated decrease of S coincides 
with the accelerated increase of jmax already observed. 
The change in i5 can be associated with a concurrent 
change in the functional form of the magnetic energy 
spectrum. After t « 2.5, the energy spectrum in the 
2048'^ run shows a slight excess of energy at small scales 
compared with Eq. p5p . Indeed, the functional form for 
the energy spectrum given in Eq. p5|l may become too 
simple after t w 2.5. We have also verified that the ac- 
celeration in the decrease of 6 persists when the range of 
wavenumbers used for the fit is changed. 

To evaluate the credibility of the computation, we set 
as our criterion for accuracy 5k„iax ^ 2, where 5{t) is fit 
over the range k £ [2, 400] for the (2048'^) run, which has 
maximum wavenumber kmax = N/3 = 682 after dealias- 
ing. By this measure, the computation is reliable up 
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to i PS 2.61, in agreement with the above examination 
of ahasing levels in the current density in configuration 
space. For each run. the threshold of reliability based on 
the Skmax > 2 criterion is marked by a dotted horizontal 
line in Fig. [2l that is, each line (labeled by its spatial 
resolution) corresponds to the minimum value of S that 
can be resolved by each simulation. Therefore, the point 
at which S{t) intersects the horizontal cutoff line for the 
given resolution marks the last reliable time. It is critical 
to note that the 2048^ run shows the acceleration in the 
evolution of 5 taking place well before the cutoff time, 
while the other runs show no indication of acceleration 
prior to their respective cutoff times. 

A third criterion for ensuring accuracy of the compu- 
tations is the conservation of the invariants of the MHD 
equations. The time at which conservation breaks down 
in the simulations is in good agreement with the time 
obtained from the two other criteria described above. To 
give an example, we note that in the 2048^ run the energy 
conservation begins to break down (at the level of 5 x 10^^ 
of the initial energy) at t ^ 2.6, which is consistent with 
the time at which S reaches the smallest resolved scale of 
the flow. 

Finally, the evolution of the spectral index n in Eq. 
as a function of time is shown in Fig. [21 for each 
run up to its respective cutoff time. For all resolutions 
we observe that the spectral index is smaller than in the 
Euler case [l^l , leading to a shallower spectrum and likely 
corresponding to sharper structures in the MHD flow. 
This could also be related to the presence of a magnetic 
field that provides a slower time scale for the spectral 
energy transfer 0, [sil, [s^ • 



D. Spatial structures in the ideal case 

We now turn to an examination of the structures that 
develop, concentrating on the magnetic field and current 
in the 2048'^ run, for which the temporal evolution shows 
an abrupt change at i « 2.5. Because of the symmetries, 
the strongest gradients tend to develop near the insulat- 
ing walls, where thin current sheets appear. The roll-up 
of these sheets, as observed in a dissipative simulation 
at high resolution of a different flow [2] does not occur 
here; presumably the reason is that the conditions for a 
Kelvin-Hclmholtz-like instability are not present, due to 
the symmetries satisfied by the vorticity and the current. 

Two strong current sheets appear near the edges of 
the box, approximately halfway up along one direction 
(z « 7''/2), whose geometry and evolution we describe 
here. Fig. [3] shows the structures in a subvolume with 
approximately 200 x 200 x 120 grid points, as well as the 
magnetic field lines, colored according to their strength. 
The current intensity is concentrated in the two thin 
sheets, which are connected to two other thin sheets 
running along the insulating walls (only one of which is 
shown in the figure). At each of the edges of the box, the 
two horizontal sheets are closest near the edge and farther 




FIG. 3: (Color online) Top: Current density intensity in a 
200 X 200 X 120 subvolume of the 2048^ run at t = 2.5. Note the 
two current sheets approaching each other. Bottom: Magnetic 
field lines in the same structure (viewed from the back). The 
current intensity in a slice is given as a reference. To the right 
and left of the slice, the magnetic field is strong (purple color), 
whereas in the transition region it decreases to ~ 1/6 of the 
maximum (yellow), corresponding to a strong local drop in 
intensity. Note also the rapid spatial rotation of the field 
lines between the two current sheets. 



apart towards the interior. The magnetic field around 
these structures is smooth and large-scale, strong outside 
the sheets but weak between the two sheets; the mag- 
netic pressure gradient thus provides a force that moves 
the sheets closer together. Each individual sheet cor- 
responds to a sheared magnetic field configuration, and 
the field between the sheets undergoes a rapid rotational 
shear. 

As time evolves, two distinct processes take place: each 
sheet becomes thinner, and the two sheets continue to 
be pushed together by the magnetic pressure; these pro- 
cesses occur at different rates. Before t w 2.5, the thick- 
ness of the sheets eg decreases faster than the separation 
distance ca between the two. The magnetic field lines 
between the sheets rotate smoothly from the direction 
observed above the top current sheet to the direction be- 
low the bottom current sheet. At t ~ 2.5, es and ca 
become comparable, and subsequently the rate of near- 
collision of the two sheets surpasses the rate at which 
each sheet thins, hence the acceleration in the evolution 



6 



of the logarithmic decrement S and in the growth of jmax ■ 
As stated previously, the computation has to be stopped 
when either es or becomes comparable to the distance 
between grid points Ax. 

The rotational shear intensifies as the sheets approach 
each other. Indeed, each magnetic field line in Fig. [3] is 
separated by two grid lengths in the cross section. At 
the time shown in the figure {t = 2.5), a sharp spatial 
rotation of the field lines of approximately 7r/2 can be 
seen to occur over a short distance, corresponding to a 
few field lines (about 8 grid lengths), while the drop in the 
amplitude of the magnetic field between the two sheets 
occurs over about twice that distance. The transition 
region becomes thinner rapidly after t = 2.5, pushed by 
the larger magnetic pressure outside the two sheets, and 
by the time the simulation is stopped the magnetic field 
rotates by 7r/2 almost abruptly. 

In comparison, structures in the vorticity are 
smoother. This is consistent with the fact that the accel- 
eration of 5 is not observed when the logarithmic decre- 
ment is computed from the kinetic energy spectrum. 



Spectro ot t = 2.50000 
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FIG. 4: (Color online) Total energy spectrum for the IMTG 
flow on a grid of 512"^ points in the ideal case; results are 
obtained using a symmetric code (diamonds) and a run with 
no imposed symmetries (solid line). No difi^erence is observed 
except, at a low energetic level, in the lowest mode of the 
computation. 



E. Comparison with a run without imposed 
symmetries 

One question to be addressed is whether the imposed 
symmetries alter in substantial ways the dynamics of the 
flow. Or, in other words, if the symmetries of the ini- 
tial conditions are preserved by the time evolution when 
symmetries are not imposed in the simulation (at least 
for times long enough to justify the savings in memory 
and computing time in free decaying simulations). It 
is known for example that imposing symmetries can slow 
the growth of the MHD dynamo instability for some flows 
fl^ , in which the fastest growing mode does not obey the 
imposed symmetries. To quantify the effect of imposing 
the symmetries in the evolution of the IMTG flow, we 
performed two simulations starting from the same ini- 
tial conditions and with a resolution of 512'^ grid points: 
one with imposed symmetries, and another without any 
imposed symmetry. 

Fig. |4] shows the energy spectrum near the end of the 
computations {t = 2.5). Diamonds correspond to data 
produced with the symmetric code, and the solid line 
corresponds to the equivalent run with a general code, 
which allows any non-symmetric mode to evolve freely. 
No difference is observed except at the lowest mode and 
at an energetic level close to round-off errors. The ampli- 
tude of any of these modes can be used to estimate how 
much energy is contained in modes that do not satisfy 
the symmetries. In particular, the energy in the modes 
in the fc = 1 shell: k = (1,0,0), (0,1,0), and (0,0,1), 
which do not satisfy the symmetries of the IMTG flow, 
in the run without imposed symmetries is 3 x 10~^^, of 
the order of the truncation error. Note that oscillations 
between even and odd wave numbers are visible in the 
energy spectra of the kinetic Taylor-Green vortex. They 



may be caused by interferences between structures lo- 
cated near the symmetry planes that are separated by 
a distance of tt. These oscillations can be removed by 
averaging the spectra on shells of width Ak = 2 [l^ . 

We also verified through visualizations that the physi- 
cal structures that appear in the flow are nearly identical 
between the two runs. Therefore, we expect no difference 
in behavior in the ideal case for a full run without im- 
posed symmetries up to the resolution and time achieved 
here and conclude that a run with imposed symmetries 
gives the same results at a substantially lower cost. 



IV. CONCLUSION AND PERSPECTIVES 

The existence of singularities of fluid equations has 
been a subject of many studies employing symmetric 
flows, since evidence of their occurrence in finite time 
is all that is needed to disprove the regularity of the 
equations, however constrained the flow in question may 
be. In this sense, use of the IMTG flow proposed here 
is a well-accepted procedure, and it allows for reaching 
resolutions never attained before in MHD using spectral 
accuracy. 

Using an IMTG initial configuration and maintaining 
the symmetries at all times, we found that the ideal flow 
appears to remain regular throughout the interval of reli- 
ability, although an acceleration of the exponential decay 
occurs just before the limit of the resolution is reached in 
the run corresponding to a regular grid of 2048'^ points. 
Our conclusion is based both on examination of the tem- 
poral evolution of current and vorticity norms and on the 
logarithmic decrement fitting technique. Higher resolu- 
tion runs will have to be performed to check later time 
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behavior — in particular to confirm whether a rotational 
discontinuity spontaneously appears. The study does not 
preclude the occurrence of singularities at later times, nor 
does it preclude their occurrence in other types of ideal 
MHD flows, as considered in pioneering studies making 
use of flnite difference adaptive mesh refinement ml o r 
specific configurations of a weakly compressible flow [29| . 
Another possibility for computing longer in a DNS might 
be to use a filter for small scales different from the stan- 
dard 2/3 deahasing rule, as proposed in p6l |. 

Apart from the singularity issue, we saw an accelera- 
tion of the dynamics never observed before in this type 
of study, either in an Euler flow or in ideal MHD. This 
phenomenon can be linked specifically to the geometric 
configuration of two current sheets becoming thinner in- 
dividually in the presence of a large-scale shear and, at 
the same time, accelerating toward each other because of 
a differential in magnetic pressure. 

The next obvious step is to compute the turbulent de- 
cay of the IMTG fiow at a magnetic Prandtl number 
of unity and examine the interactions between turbu- 
lent eddies and Alfven waves 0; of particular interest 
is the potential emergence of non-Kolmogorovian energy 
spectra. For example, it was predicted on a phenomeno- 
logical basis by Iroshnikov and Kraichnan that Alfven 
waves would slow the dynamics of MHD flows, relative 
to neutral fluids, and thereby produce a shallower spec- 
trum; this shallower spectrum can also be derived using 
the Lagrangian renormalized approximation generalized 
to MHD [3l|. Other spectral laws have been derived, for 



instance, in the context of anisotropic ("weak") MHD 
turbulence; scaling laws depending on the characteris- 
tic timescales of the system may emerge as well, such as 
when forcing functions are introduced with their intrin- 
sic correlation times in the presence of a uniform strong 
magnetic field [33|. By using the IMTG fiow and imple- 
menting numerically its symmetries to be able to attain 
substantially higher numerical resolutions, we can con- 
tinue the search for scaling laws in MHD. 

Finally, the properties of the second TO fiow (corre- 
sponding to the conducting case) will also be examined in 
a future work. Since rotational structures, including dou- 
ble discontinuities [1^, have been observed in interplan- 
etary space in a much more complex environment involv- 
ing compressibility and Hall currents, the fiows proposed 
here may prove useful in understanding the interplay be- 
tween turbulent eddies, waves, and small-scale structures 
in interplanetary space and the heliosphere. 
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University of Colorado, and an IBM Shared University 
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part from an NSF ICERT Fellowship ( Joint Program in 
Applied Mathematics and Earth and Environmental Sci- 
ence) at Columbia. P.D.M. is a member of the Carrera 
del Investigador Cientifico of CONICET. Finally, visual- 
izations use the VAPOR software for analysis of terascale 
datasets fsJ l. 
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